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We introduce an effective quark-meson-nucleon model for the QCD phase tran¬ 
sitions at finite baryon density. The nucleon and the quark degrees of freedom are 
described within a unihed framework of a chiral linear sigma model. The decon- 
hnement transition is modeled through a simple modification of the distribution 
functions of nucleons and quarks, where an additional auxiliary held, the bag held, 
is introduced. The bag held plays a key role in converting between the nucleon and 
the quark degrees of freedom. The model predicts that the chiral and the deconhne- 
ment phase transitions are always separated. Depending on the model parameters, 
the chiral transition occurs in the baryon density range of (1.5 — 15.5)no, while the 
deconhnement transition occurs above 5no, where no is the saturation density. 


PACS numbers: 12.39.Fe,12.38.Mh,25.75.Nq 
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INTRODUCTION 


The basic problem of QCD thermodynamics is to understand the conversion from 
hadrons to quarks and gluons and how this is related to the underlying chiral and 
deconhnement transitions. This problem has major relevance to the physics of heavy 
ion collisions and compact stars At present there are several proposals to effec¬ 

tively convert between the hadronic and the quark-gluon phase, with partial success. 
At hnite T, quark degrees of freedom can be suppressed by the Polyakov loop Bfli. 
Mechanisms for excluding composite degrees of freedom are the spectral function 
method js-^ and the excluded volume method, for recent works see Refs. 14- 
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In this work we construct an effective quark-meson-nucleon model with two flavors 
for the QCD phase transitions at hnite density. We place several restrictions to our 
approach: hrst we must take into account both the nuclear and the quark degrees 
of freedom. Second, the model should respect global symmetries of QCD, i. e. 
the chiral symmetry and the scale invariance. Third, quark degrees of freedom 
must be excluded at the nuclear matter density, and nuclear degrees of freedom 
must be excluded at some high density. We require that the model reproduces 
the nuclear matter ground state. Finally, a unihed description must encompass 
couplings of nucleons and quarks to the same bosonic mean-helds generated from a 
unique vacuum potential. 

For nucleons, both chiral and scale invariance can be accommodated in the parity 
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model 
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24| 25|-[27| . The quark sector is described by a linear sigma model 


29[ | coupled to the dilaton 


30|. We hnd that the small a mass as required for 


reasonable nuclear matter properties 2l(] results in a shallow potential in the a 
direction. As an immediate consequence, pure quark matter appears at baryon 
chemical potentials below the value of the vacuum nuclear mass. In order to 
solve this problem we generalize the idea of statistical conhnement from effective 
model studies at hnite T to hnite densities. While at hnite T the Polyakov loop is 
used to statistically suppress thermal quark huctuations j^, its extension to hnite 
baryon chemical potential is problematic. In this work we use instead the concept 
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of infrared confinement {m, 321 in order to modify the Fermi-Dirac distributions 
of quarks. We consider a simple model where the Fermi distribution of quarks is 
restricted to momenta above b, where 6 is a new auxiliary held in our model, which 
we name the bag held. The hnite value of b in the vacuum and at low temperatures 
and densities is guaranteed by a new phenomenological vacuum potential. We ht 
the parameters of this potential to the QCD vacuum energy and by matching the 
pseudo-critical temperature for the chiral and deconhnement transition known from 
the lattice simulations at Nf = 2 33|, [3^. The model is compared to a Polyakov- 
quark-meson model by pointing out their similarities and diherences. 

Next, we introduce a generalization of the distribution functions for the nucleons 
in such a way that their Fermi surface is restricted only to low momenta 32|, i. e. 
below some value ab, where a is a new parameter. This leads to a construction 
of a combined quark-meson-nucleon model. We hnd that the minimization of the 
thermodynamic potential in the b held acts to convert the nucleons to quarks as the 
baryon chemical potential is increased. One of the main consequences of this model 
is that the chiral and the deconhnement transitions at T = 0 are separated. 

This paper is organized as follows: in Sec. [TTl we introduce the nucleonic model 
and briehy consider the nuclear matter ground state. Sec. Hill is devoted to the quark 
degrees of freedom. Here we introduce the 6-held potential and discuss statistical 
conhnement of quarks. The Sec. HV] describes a combined quark-meson-nucleon 
model. The main results of this paper are given in Sec. |Vl while in the following 
Sec. I VI I we make our conclusions. In the Appendix we solve a simplihed version of 
the quark-meson-nucleon model. 


II. NUCLEONIC SECTOR: PARITY DOUBLET MODEL WITH 

DILATON 


We consider the Nf = 2 nuclear parity doublet model within the mirror assign¬ 


ment of chiral symmetry 


18l-l24l| 2^^. We prefer to use a linear realization of 


the chiral symmetry which allows the description of chiral symmetry restoration. 
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The restoration of chiral symmetry in QCD dictates that hadrons of opposite parity 
become degenerate, but not necessarily massless. A hnite chirally invariant mass is 
then modeled by a parity doublet model with mirror assignment. Note that the pos¬ 
sibility of a chirally invariant contribution to the nucleon mass was recently hinted 


by lattice calculations 
at hnite T 


35| . First lattice QCD simulations of nucleon parity partners 


36l | hnd that their masses become degenerate by crossing the QCD phase 


transition. 

The nucleonic part of the model Lagrangian 


Q,Ei 


coupled to the dilaton 


25 


27| 


IS 


Cn = + C/xX(^175^2 - ^275^1) 

-h -h iysT ■ 7 r)Ti g2^2{o- — i'y^T ■ 77)4/2 

- - g^^20^2 ■ 


( 1 ) 


where are the nuclear chiral partners. The fermions 4/i^2 are coupled to the 
chiral helds (cr, tt), to the held and to the dilaton y;. The mass eigenstates are 
given as 



1 / 6-5/2 ^^g-5/2 


where 


with masses 


\/2 cosh 5 y^ge-'5/2 —£<5/2 

gi + g2 0- 



( 2 ) 


sinh5 = 


2£/x X ’ 


mN± = 


2 L 


( 3 ) 


(^1 + ^2) V2 + ^ 

The state N+ is the nucleon N{938) while is its parity partner conventionally 
identihed with A^(1500). The meson contribution is as follows 


Cm = 

- K - 14 - 14 , 


( 4 ) 


where 
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( 5 ) 
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T/ 2 2 

=- 


( 6 ) 


and 


^ 4 


A 

Xo 


log 


(i)'- 


(7) 


The total Lagrangian Cjy + Cm is chiral and scale invariant. All the masses in the 
model are generated by the condensation of the dilaton held in the vacnnm xo- We 
can hx by = A^Xg, where = 783 MeV. The parameters A, A^ and e are 

related to the sigma and pion masses and the pion decay constant f-,, as 


A = 


m: 


m: 


Ax = 


m: 


3m^ 


e = 


Xo 


( 8 ) 


2/.^ ’ 2xg 

with = 138 MeV and = 93 MeV. We take m+ = 938 MeV, m_ = 1500 MeV 
21|. The nnclear matter gronnd state can be obtained by hxing the parameters 
iTT'a, 91 , 92 , 9 x, 9 uj- The dilaton potential 1^ is hxed by identifying the lowest 
glneball mass with the dilaton mass = 1700 MeV |37|, ISSj, and by hxing the 
valne of the glnon condensate. The conventional valne of the glnon condensate 
~ (331 MeV)"^ 39| is accompanied b y la rge nncertainties, {^G^yG^'') ~ 
(300 — 600 MeV)"^. For recent acconnts see 40l-l42| and references therein. By the 
relation for the trace anomaly this can be translated into the following range for the 
QCD vacnnm energy Cvac — (193 — 386 MeV)"^. Assnming that the QCD vacnnm 
energy is dominated by the dilaton potential, the glnonic bag constant is estimated 
to 5 ~ (273 — 546 MeV)^. From the relation for the dilaton mass 


2 


AB 

^ ’ 
Xo 


(9) 


we obtain xo — 87.79 — 351.17 MeV. 

Since the dilaton is heavy = 1700 MeV, it practically does not inhnence 
the nnclear gronnd state and we can adopt the model parameters from 2l|. This 
hxes rria = 370.63 MeV, 91 = 13.00, 92 = 6.97, 9 y. = 4.39 and 9 ^^ = 6.79. The 
corresponding thermodynamic potential in the mean-held approximation is 


= 14 + lA + W 


X=N± 


Q 


X 


( 10 ) 
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= 7x 


d^p 

(27r)3 


[Tlog(l-/x) + Tlog(l-/x)] , 


( 11 ) 


where the functions fx are the Fermi-Dirac distributions 

1 ^ 1 


fx = 


fx = 


and Ex = a/p^ + px = I^b — Qlj^i with 'yx = “2 x 2 being the spin-isospin 
degeneracy factor. We minimize the potential with respect to a, x and uj 


d^l X 2 X 

— = a +An 




E dmx „ 
- 0 > 


X=N± 


da 


dfl ^ 2 

^ = -Aa;X Uj + guj 2^ Px = 0 , 

^ X=N± 

X 2 |'^X3 o \ 2 |dA^i 

_ = -A^crx+—X -2eax-A^xa; -^5—log 
OX ^ Xo 

E dmx 

= 0 , 


( 


Vxo 


( 12 ) 

(13) 

(14) 


X=N± 


dx 


where the scalar and the baryon number densities are, respectively 

dfp mx 


and 


sx = In 


Px = In 


(2x)3 Ex 
f d?p 


{fx + fx) , 


(2x) 


{fx - fx) ■ 


(15) 


(16) 


As mentioned above, with the present parametrization, the contribution of the 
dilaton held to the ground state properties is numerically negligible. One can check 
whether y can inhuence the nuclear matter equation of state at all. In principle. 


this is possible but with a lower dilaton mass, see 


^|. With a mass of = 1700 


MeV used in this work the impact of a dilaton held, e. g. on chiral restoration, is 
expected only at much higher densities. We briehy discuss this possibility in Sec. |Vl 
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FIG. 1: Quark number density in the QM and in the bQM model. 


III. QUARK SECTOR: LINEAR SIGMA MODEL AND STATISTICAL 

CONFINEMENT 


We introduce the quark-meson (QM) coupling according to the linear sigma 
model 

Cq = iq^q + gqq{a + i'y^T ■ 7r)q - W , (17) 

with the same parameters in the cx-potential as in Eq. ([5]) where we considered 
nuclear matter. For the quark-meson coupling we use Qq = 300 MeV//^. The 
dilaton is too heavy and hence not essential for the following discussion, so we set 


X = Xo by hand. The thermodynamics of this model was studied in Refs. 


44 


46| 


A. Shallow potential and early onset of quarks 


Using the same model parameters (i. e. low rria) that are constrained in Sec. [TTl 
as to obtain reasonable nuclear matter properties gives a shallow potential in the a- 
direction. A shallow potential leads to a chiral phase transition in cold quark matter 


at rather low densities, even below uq 


45[ |. The corresponding equation of state has 


a zero-pressure point where quark matter has a hnite density. This is illustrated by 
the red line in Fig. [1] where we plot the density of quarks pq as a function of the 
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quark chemical potential Hq = /is/S. In our present model we obtain 3fiq ~ 750 
MeV which is lower than the nucleon mass. 

This is a striking problem: quarks appear too early due to a shallow potential. A 
common way to solve this problem is to adjust the chiral potential in the quark sector 
independently of the nucleonic sector. Such a treatment is not possible in a unified 
description of nucleon and quark matter. Also, note that a larger m„ is not favored 
in view of nuclear matter ground state. For example, with typical values in quark 
models ^ 600 MeV, compressibility of nuclear matter at saturation increases by 


f-Q 


the 


an order of magnitude from its experimentally suggested range j2]|. In Ref. 
problem of a flat potential was circumvented by assigning a bare mass term of 200 
MeV to the quarks. We conclude that the essential missing physics is confinement 
of quarks in the infrared region where their interaction becomes strong. To avoid 
this inconsistency below we introduce a simple model of statistical confinement. 


B. Statistical confinement of quarks 

The concept of statistical conhnement is very successful at finite temperatures 
where the Polyakov loop is used to modify the quark distribution functions j^. At 
finite density and small temperature the center symmetry is badly broken so we 
cannot use the Polyakov loop. As an alternative, we propose a modification of the 
quark distribution functions via the following ansatz 

Uq = - b‘^)fq , Hq = - b‘^) fq , (18) 

where 6 is a parameter, and fg and fq are 

f = _ I _ f = _ I _ 

the Fermi-Dirac distribution functions for quarks and antiquarks, respectively. Ob¬ 
viously, in Eq. (1T8|1 quarks with momenta p^ < 6^ are suppressed. 

This is one possible way to restrict thermal quark fluctuations at low momenta. 
It is similar to the concept of infrared confinement used in the Dyson-Schwinger 
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vacuum studies j^, ^ and in the NJL model ^ l49l-l5l|. The infrared cutoff 


is in-line with the idea of in-hadron condensates 
the holographic hard wall 


52| . and is also implemented in 


53| and soft wall ^ models. Intuitively, 1/6 can be 


understood as a typical size of a hadron, so that due to the uncertainty principle 
quarks cannot have momenta lower than b. 

With a sharp cutoff in the distribution function it is not possible to saturate the 
Stefan-Boltzmann limit at high temperature and/or density. Essentially, b must be 
a medium dependent quantity. A thermodynamically consistent way to achieve this 
is to promote 6 to a field generated by some potential 14- The minimization of the 
thermodynamic potential in the 6-direction results in 6 being a medium dependent 
quantity. Since the potential 14 is an additional contribution to the bag pressure, 
this prescription can be understood as a self-consistent way to generate a medium- 
dependent bag pressure. We therefore name this model the bag-quark-meson model 
(bQM), and the field 6 is named the bag field. We consider the bag field as a non- 
dynamical, auxiliary held. This held is responsible for statistical conhnement, in 
spirit similar to the Polyakov loop. 

Taking into account the modihcation of the distribution functions, and the vac¬ 
uum potential 14 ; the thermodynamic potential of the model becomes 




where 


~ Iq 


d?p 


[riog(l-n,)+ riog(l, 


(19) 


( 20 ) 


with Hq given by (ITSll and 7 ^ = 2 x Wj x W = 12 for two havors. The gap equations 
are 


do, 

= -Kxlo- + - exl + QqSq = 0 , 


and 


where 


dXl dVb 


db db 


ZUq 0 , 


( 21 ) 


( 22 ) 


d^p 


rUn 


(27r)3 Eq 


{Uq -F Uq) . 


Sq=lq 


(23) 
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FIG. 2: The mean-fields a and h in the bQM model as a function of T at /ig = 0. Left 
panel: results for different values of the parameter Kb, with Cvac = (250 MeV)^. Right 
panel: results for Cvac = (250 MeV)^, Kb = 155 MeV and Cvac = (386 MeV)^, Kb = 400 
MeV. In this case, the pseudo-critical temperature of the chiral transition is fixed to 
Tc = 170 MeV. 

The term zUg is a boundary contribution of (I2UD 

= Iq^ [^10g(l - fq) + Tlog(l - /g)] p 2 =fc 2 • (24) 

The goal is to suppress quarks at low T and/or pg with a large value for h, and 
to have lower values of b at high T and/or pg. Remarkably, this is accomplished 
through the minimization of the thermodynamic potential with respect to b: the 
thermal correction zUg to the gap equation for b acts to reduce b. 

C. Parametrization of the bQM model 

Unlike the Polyakov loop potential, the potential Vb cannot be constrained by 
symmetry, so one should try different forms. Below we consider a special case with 
a minimal number of terms that yields a hnite value of b in the vacuum. Namely, 
we choose 

14 = ^ _ 


(25) 
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where the new parameters Kb and Xb must be determined. The non-trivial vacuum 
expectation value of this potential is Bq = ^JK^/Xb. One can further motivate this 
choice as follows. Since the b held generates statistical conhnement in our model, 
we associate it with the chromo-electric part of the gluon sector. The essential 
reason for this is that hnite T lattice computations show that the chromo-electric 
part of the gluon condensate drops across the pseudo-critical temperature |55l457|. 


Also, the chromo-electric sector is governed by the zeroth component of the gluon 
helds, like the Polyakov loop or the quark-antiquark potential, so it is essentially 
non-dynamical. In that sense, such an identihcation can be considered natural. 

Since the 6-potential fl25|) yields an additional contribution to the vacuum energy, 
we must ensure that the total vacuum energy remains correctly saturated. We will 
consider a particular case where the vacuum energy created by the y-held and the 
6-£eld are equal in magnitude = 14p, so that a half of the total vacuum energy 
has its origin in the chromo-electric sector modeled here by I 4 . In this 

language, the dilaton y natu rally represents the chromo-magnetic component of the 


gluon condensate 


mm 





59l |. which survives the chiral transition. We thus take 
4 


Kl 


T/ _ ■''b 


^vac 


(26) 


where the total vacuum energy is Cvac — (193 — 386 MeV)"^. We hxed a value for 
Vbg (within the range for Cvac) and then changed Kb. We found that in general low 
values of Cvac and low values of Kb produce a hrst order chiral transition in the bQM 
model a.t fig = 0 and hnite T, see the left panel Fig. [2] for Kb = 50 MeV. Since 
lattice QCD shows that the chiral transition is a crossover this sets a lower bound 
on Kb- Increasing Kb while holding Cvac fixed leads to a decrease of bo. The result is 
that the chiral transition turns into a crossover with its onset shifted towards lower 
temperatures as shown on the left panel of Fig. O 

As a second constraint, we choose the pseudo-critical temperature as ob¬ 
tained in lattice QCD simulations. For Nf = 2 -|- 1 with physical quark masses 
the Wuppertal-Budapest collaboration obtained Tr = 147(2) (3) MeV [6^, while the 
HotQCD collaboration quoted = 154(9) MeV [Mi- Since in this paper we work 
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33, 3- 


with Nf = 2 we will use a slightly larger value Tc ~ 170 MeV [33|, l34| . This gives 
^vac = (250 MeV)^ and Kf, = 155 MeV. From (1261) we hud A;, = 0.074. The resulting 
a and b mean-helds obtained by solving fl^ and fl22|) at hnite T and /ig = 0 are 
shown on Fig. [21 where bo = 570.3 MeV. 

With the appropriate modihcation of the quark number density 


d^q f d?P ^ ^ 

we hnd that the onset of quarks for eyac = (250 MeV)^ in the bQM model at T = 0 
is at fiq ~ 450 MeV. This numerical value must be contrasted to the one in the QM 
model, where fiq ~ 250 MeV. 

It must be stressed that the obtained parameters Kb and \b are not unique. The 
uncertainty in the QCD vacuum energy provides a range for the parameters Kb 
and Afc. The hxed value of the vacuum energy Cvac = (250 MeV)"^ is the smallest 
value which still satishes the above-mentioned constraints. Taking higher values of 
evac, while keeping Tc hxed we hnd that both Kb and \b increase. Interestingly, the 
vacuum expectation value bo = K^/Xb shows a reduction of less than 10% in the 
range evac = (250 —386MeV)^. The sensitivity to Cvac is rehected in the temperature 
at which the 6-held experiences a rapid change. Higher values of Cvac represent a 
larger energy barrier for quark huctuations so a delayed and more gradual change in 
the fe-held is expected. This is visible on the right panel of Fig. |2] where the b held 
is shown for Cvac = (250 MeV)^ and Cvac = (386 MeV)^. 


D. Similarities and differences between the bQM and the PQM model 

Since a hnite b held reduces the strength of the quark thermal huctuations, it 
is intuitively clear that it will act to increase the critical temperature associated 
with the chiral transition (see also left panel of Fig. [2]). In the chiral limit we can 
show this relation analytically by hnding a zero of the cr^ coefficient in the Landau 
expansion of the thermodynamic potential flT^ in powers of cx^. The result can be 







13 



FIG. 3: A comparison between the normalized mean fields in the bQM and the PQM 


models as functions of T at fig = 0. The bQM results are presented for evac = (250 MeV)^. 


cast in the following parametric form 


T2 



0 , 


(28) 


where 


J'(x) = ^(-3x^ + vr^ + 6xlog(l + e^)) + ^Li 2 (-e"') , (29) 

and Li 2 (x) is the polylogarithm of order 2. In the limit h/T 1 this can be 
simplihed to 


r,»Q“ ~ (PhnI + !LA . (30) 

V lq9 ) 

The second term under the square root is a correction to the usual QM model result 
for Tc- It is instructive to compare this result to the one obtained in the Polyakov- 
quark-meson (PQM) model 


rpPQM ^ f I 

" V lq9‘^ J 


1/2 


(31) 


where (j) is the background gauge held related to the Polyakov loop <I> through 
$ = [1 + 2cos(0/T)]/3, see e. g. 5]. 

In the PQM model the Polyakov loop potential saturates the transverse (phys¬ 
ical) gluon contribution to the thermodynamics. On the other hand, in the bQM 
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model the potential H does not contain the transverse glnons. This is the main 
difference between the PQM and the bQM models. In the PQM model, deconfine¬ 
ment transition can be characterized by the Polyakov loop. Then, one can dehne 
the psendo-critical temperatnres governed by the peaks of the qnark condensate and 
the Polyakov loop and interpret their approximate coincidence as obtained in lattice 
QCD simnlations 6^. However, the bQM model in its present form cannot predict 
the deconfinement transition temperatnre, nor the coincidence of the two respective 
pseudo-critical temperatures, because the b held is not a conhnement-deconhnement 
order parameter. 

To compare results for the bQM and the PQM models, we have calculated the 
a held and the Polyakov loop in a PQM model. A concrete parametrization of the 


Polyakov loop potential is taken from 


62| , while for the quark sector the same linear 


sigma model flTTl) is used. On Fig. [3] we compare our results where Cvac in the bQM 
model was set to Cvac = (250 MeV)^. We hnd that the a held in both models are 
closely matched, with peaks in da/dT at 163 MeV and 170 MeV (the htted value), 
respectively. Obviously, a more realistic consideration should include explicitly the 
bosonic excitations. 


IV. A COMBINED QUARK-MESON-NUCLEON MODEL 

In this section we construct a combined quark-meson-nucleon (QMN) model. Our 
guiding requirement is to exclude quarks at low density, and to exclude nucleons at 
high density. In the previous Section we have introduced the concept of statistical 
conhnement of quarks, and now we propose a similar but opposite modihcation of 
the distribution functions for nucleons 

nN± = - pQ/jVi , nN± = - pQ/jv^ (32) 


With this ansatz nucleon huctuations with momenta p^ > {aty are suppressed, 
where a is a new parameter of the model. A similar form has been used in 3^. In 
Section flV Bl we will determine the possible range for the parameter a. 
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A. Model setup 


The thermodynamic potential of the QMN model is obtained as the following 
sum 



(33) 


X.iV±,j 


with dehned as in flTT]) with the appropriate modihcation of the distribution 
functions fN± according to Eq. ([32]), while flq is dehned in fl20D . 

The gap equations are obtained by minimizing the thermodynamical potential 
(133]) . The u and x gap equations remain unchanged by the inclusion of quarks, being 
given by Eqs. ffT5]l and (ITT|) . respectively, and the proper replacement t n 7 v± 

according to Eq. fl32|) . Thee gap equations for a and b are modihed as follows 



(34) 


X=N±,q 



(35) 


X=N± 


where 



(36) 


and zUq is given in (12T|) . 

The pressure is calculated by evaluating the thermodynamic potential at its min¬ 
imum p = —normalized with the constant flo in such a way that the physical 
vacuum has zero pressure. The total baryon number density is 


dQ 


1 

PN+ + PN- + -Pq , 


Pb — -w— 
OPB 


where we used Pg = /is/S. The nucleon and quark particle fractions are dehned as 



(37) 


For the transition from nucleon to quark degrees of freedom the gap equation 
for the 6-held, Eq. fl35ll . is crucial. Note that, as a consequence of the Leibniz rule, 
the nucleon and quark contributions in (j35ll have opposite signs. In the low density 
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phase nucleons will favor finite b and as a consequence quarks are suppressed. The 
appearance of quarks at high densities acts to reduce b and therefore exclude nucle¬ 
ons. Therefore, Eq. fl3^ controls the relative abundance of nucleons and quarks. In 
order to better illustrate this point we have solved a simplihed model with only the 
b held in the Appendix. 


B. Parametrization 


ruo- (MeV) 

(MeV) 

B (MeV)^ 

91 

92 

9x 

9i^ 

9q 

A 

Ax 

Kb (MeV) 

Aft 

370.63 

1700 

297.30 

13.0 

6.97 

4.39 

6.79 

3.22 

6.84 

3.71 

155 

0.074 


TABLE I: The parameters of the quark-meson-nucleon model. 


The parameters of the vacuum potential (as well as the effective masses mx and 
chemical potentials /ix) are dehned in Sec. [TTl and Sec. IIIIl and collected in Table 
[B In particular, we will use 14 with Cvac = (250 MeV)"^ as discussed in Sec. IIII Cl 
The remaining parameter a is chosen so that the effective UV cutoff ab^ for the 
nucleon distribution functions does not spoil the nuclear matter ground state. This 
sets the minimal value for ab^ > 300 MeV. Assume now that at some /i^ there is 
a transition to quark matter. From the following consideration this will lead to a 
useful estimate on the upper bound on a. Consider Eq. (|35|1 in the limit of large 
/ifi. At T —)■ 0 the boundary terms read 


[cib')^ 

'^N± —JN ^^2 “ Ex±)9{^n — Ex±) , 

KJ, -*■ - E,)e(M, - E,) . 


(38) 


(39) 


We assume a, a; —)■ 0 in this limit. In addition, we consider a case where /is is large 
enough that we can also ignore the chirally invariant nucleon mass. Using 
Eq. fl2^ simplihes to 


Oil 2 t ^ ,3 

+ ^bb + 


27r2 


^ —2a^7Ar -|- ^ j /iB — 0 . 


(40) 
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(MeV) (MeV) 


FIG. 4: The mean fields for the QMN model for two different values of the a parameter 
indicated in the figure. As fiB is increased there is a clear imprint of the liquid-gas and 
the chiral phase transitions on the a field shown by the dash-dot black line. The drop in 
the b field at high /xb is associated with the deconfinement transition. 

Since 'Jn = 7q/3, the vanishing bracket defines Omax = 2“^/^. For bo = 570.3 MeV 
we have ctmax^o = 452.6 MeV. In the following results we consider two values of the 
a parameter: abo = 300 MeV and abo = 440 MeV, which is close to the benchmark 
value Omax^o = 452.6 MeV. We also calculate the transition points for several a’s in 
between. 


V. RESULTS 


Without the modification of the nuclear distribution functions the parity doublet 
model with dilaton introduced in Sec. m has a weak first order chiral transition at 
/iB — 2110 MeV for the parameters used here. This should be contrasted to the 
value of /iB — 1725 MeV found for the same parameter set j^, 22] but using a model 
without the dilaton. Such a delay in the chiral transition takes place because the 
dilaton field is essentially coupled to the density through the oj field, see Eq. ([6]). 
In the case where the distribution functions are modified we must deal with the 
uncertainty in the parameter a and Cvac- We show explicit calculations for the case 
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(MeV) |ig (MeV) 

FIG. 5; Quark and nucleon particle fractions for the QMN model for two different values 
of the a parameter as indicated in the figure. 

Gac = (250 MeV)"^ and for two limiting values of a. Finally, we will discuss sensitivity 
to the uncertainty in Cvac- 

With abo = 300 MeV the onset of chiral-symmetric phase occurs at a lower 
Hb- This can be seen on Fig. 0] showing the solutions of the gap equations in the 
QMN model. This is easily understood: in this model it is the nucleons that restore 
the chiral symmetry - quarks appear only at higher see Fig. [TJ The cutoff in 
the nucleon momenta limits their density, and the u held reaches a plateau after 
Hb — 1000 MeV. Therefore, the shift of /is due to the u held is diminished and, 
as a consequence, the chiral transition happens at lower than in the model 
without a cutoh. At aho = 300 MeV the chiral transition is hrst order and occurs at 
fiB = 1135 MeV. After this point, the a held drops almost to zero and the parity- 
doublet partners become degenerate and equally populated. In the case ah^ = 440 
MeV, the chiral phase transition is a crossover. From the peak in da/dfiB we extract 
fiB = 1473 MeV. At that point the mass splitting of the parity partners is around 
18% of its chirally invariant contribution in the vacuum given by g^Xo = 790 MeV. 

At some higher /i^ the 6-held decreases, as can be seen on Fig. 01 The reduction 
of the 6-held suppresses nucleons simultaneously enhancing quark huctuations, ac¬ 
cording to Eqs. (I3^ and (081) . respectively. Therefore, this marks the deconhnement 
transition in the QMN model. This is demonstrated on Fig. [5] where we plot the 
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(MeV) 

FIG. 6: The baryon number density ub (in units of saturation density no) as a function 
of the baryon chemical potential /is in the QMN model and for two different values of the 
a parameter indicated in the figure. 


nucleon and quark particle fractions defined in Eq. flHTll . 

The strength of the transitions strongly depends on the value of a. For abo = 
300 MeV the deconhnement transition is accompanied by a jump in the b held - 
consequently the baryon density has a jump, so it is rightful to consider this as 
a hrst order phase transition. On the other hand, for abo = 440 MeV, the 6-held 
reduces gradually so that the deconhnement transition is in fact a crossover, allowing 
for a wide region of a mixed phase of nucleons and quarks, see Fig. [5l In this case 
we conventionally mark the point of the deconhnement transition with at which 
Yq = Yn± holds. The model predicts that the deconhnement transition happens at 
/i^ = 1716 MeV for the case abo = 300 MeV, while it occurs at /i^ = 2130 MeV, 
for abQ = 440 MeV. The diherence between these two cases is seen for example in 
the behavior of particle fractions on Fig. [5l 

The chiral and the deconhnement transition are rehected in the behavior of the 
density ub as a function of fiB, plotted on Fig. [6], and the equation of state, plotted 
on Fig. [7| in the p — ub plane. While for both values of a the chiral transition is 
visible as a small density change, for abo = 300 MeV the deconhnement transition 
has a pronounced density jump in the region hb — Q — lOno- On the other hand. 
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FIG. 7: The equation of state for the QMN model for two different values of the a. 
parameter indicated in the figure. The locations of the chiral and the deconfinement 
transition are approximately marked. In the case of abo = 440 MeV, the arrow marks the 
point where + Y]y_ = Yq. 


with abo = 440 MeV the deconhnement transition is continnons and there is no 
clear imprint on the eqnation of state. 

In Table im we collect the nnmerical valnes for the chemical potential and the total 
baryon density at the onset of the chiral and the deconhnement phase transition 
for several valnes of a. With increasing a, both the chiral and the deconhnement 
transition are shifted to higher fis- The onset of the chiral phase transition is limited 
by the valne fiB — 2100 MeV, which is the resnlt for any a > Omax- On the other 
hand, the deconhnement phase transition will be pnshed to ps —)■ cxd as a ~ ctmax is 
approached. Beyond ctmax there is no deconhnement transition. To snmmarize, we 
have fonnd that the chiral and the deconhnement transitions in the model do not 
coincide for any choice of a. 

Lattice QCD results predict that the chiral and the conhnement transitions ap¬ 
proximately coincide along the hnite T axes of the QCD phase diagram. We have 
performed calcnlations at two non-zero temperatnres, T = 50 MeV and T = 100 
MeV to test whether snch a tendency may be observed within the QMN model. At 
T = 0 MeV the separation between the two transitions is ronghly A/i^ ~ 580 — 660 




21 


abo (MeV) 

pf (MeV) 

nf (no) 

order 

(MeV) 

n% (no) 

order 

300 

1132 

1.7 

1st order 

1716 

5.8 

1st order 

350 

1220 

2.8 

1st order 

1851 

9.0 

1st order 

400 

1348 

4.4 

crossover 

1931 

11.7 

1st order 

440 

1473 

6.1 

crossover 

2130 

18.1 

crossover 


TABLE II: The baryon chemical potential and the baryon number density at the 
onset of chiral and the deconfinement phase transitions for several values of a in the QMN 
model. We also denote the order of the chiral and the deconfinement transition for both 
cases. In the case of a first order transition, the transition is defined as a lower value of 
the number density jump. In the case of a crossover, is defined by the peak in da/d^B, 
while is defined as the point where I 7 V+ + Yn- = Yq- The values of are given in 
units of the saturation density no = 0.16 fm“^. 

MeV (see Table IH), within the explored a range. With T = 50 MeV we find 
A/is ~ 560 — 880 while at T = 100 MeV we hnd A/r^ ~ 420 — 850 MeV, i. e. there 
is no clear trend in the critical parameter. 

Taking into account the uncertainty in Cvac = (250 — 386 MeV)^ does not lead to 
a significant modification in the above results. We have explicitly checked that for 
Cvac = (386 MeV)^ the b field changes more gradually with increasing (similar 
to the hnite T results presented in Fig. |2]) thereby softening the chiral and the 
deconfinement transition. 


VI. CONCLUSIONS 


The properties of matter at large baryon densities are almost exclusively consid¬ 


ered either in purely nucleonic mode 
attempts to a unihed description 


s, or within purely quark models, with rare 
32 I. I63l-l65|. A crucial ingredient for a uni¬ 
fied description should be an effective mechanism to exclude quarks in the dilute 
hadronic matter, and to exclude nucleons at asymptotically high densities. While 
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the former has been accomplished by coupling quarks to the Polyakov loop, it is 
applicable only at hnite temperatures and low densities. In this work we have used 
a simplistic modihcation of the distribution functions of both the nucleons and the 
quarks that provides a mechanism to exclude quarks at low density and nucleons at 
high density. 

We have considered nuclear matter in a parity doublet model with mirror as¬ 
signment coupled to the dilaton held. We have argued that the chiral potential 
which is htted to the nuclear ground state properties is shallow (i. e. the a meson is 
light). When the same potential is used for quark matter it yields an early onset of 
quarks. This is especially acute at zero temperature, where we hud that with such 
a chiral potential the onset of quarks happens at baryon chemical potential below 


the nuc 
66 


eon mass. Chiral models of nuclear matter usually favor low m^, see e. g. 
681], so the problem of the shallow potential seems to be somewhat general. 


However, we stress that it is a problem only if we choose to couple the quarks to 
the same a held. 

We have proposed a possible solution of this problem by introducing the concept 
of statistical conhnement for quarks. The quark distribution functions were modihed 
in such a way that the quarks are suppressed below some particular momentum. 
We associate this minimum momentum with an auxiliary scalar held, named the 
bag held. This held brings a new contribution to the total vacuum potential. We 
have found that the minimization of the thermodynamic potential in the b direction 
provides a thermodynamically consistent framework. Moreover, the b held is hnite 
in the vacuum and is reduced as temperature or density is increased. At hnite 
temperatures the mechanism of suppression of quarks is similar to the ehect of the 
Polyakov loop. We have proposed a phenomenological form for the bag potential 
and htted its parameters to the QCD vacuum energy Cvac and to Tc at = 0 known 
from the lattice simulations. Within this scheme it was possible to solve the problem 
of a shallow potential. 

Further on, we have generalized the nucleonic distribution functions, by restrict¬ 
ing their ultraviolet momentum space up to a value ab where a is an additional 
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parameter in the model. We have used the same bag potential to construct a com¬ 
bined quark-meson-nucleon model. The vital feature of this model is that both the 
quarks and the nucleons are coupled to the same bosonic helds. Then the fate of 
the QCD phase transitions was investigated at hnite densities. 

We explored a range of values for the parameter a and found that an increase in 
a delays the onset of the chiral and the deconhnement transitions. The requirements 
that the nuclear ground state is not affected by quarks, and that the deconhnement 
transition happens at some density, restrict a to a hnite range. Taking into account 
the uncertainty in e^ac within this range of a we hnd that the chiral transition 
occurs at ~ 1100 — 2100 MeV. For the deconhnement transition we predict a 
lower bound > 1700 MeV. In terms of the density the corresponding values are 
~ 1.5 — 15.5no and > Srig. For the lowest value Cvac = (250 MeV)"^ order 
of both phase transitions depends on the value of particular a. With low values of 
a both transitions are hrst order, while for the higher values of a they are smooth 
crossovers. Especially, the deconhnement transition proceeds in a broad crossover 
mixed phase where nucleons and quarks coexist. These features persist as we vary 

^vac‘ 


In this model the chiral and the deconhnement transition are always separated. 
The chiral transition is driven by the nucleonic huctuations for any a. It is suggestive 
to consider this result in the light of the calculation of the Wilson and the Polyakov 
loop on the lattice with Dirac zero modes artihcially removed. In Refs. 69|, [TOj it was 
found that in this chirally “unbroken” phase the Wilson loop still displays an area 
law, and that the Polyakov loop is almost zero. The separation of the transitions 


might be considered as a manifestation of the quarkyonic phase 


71 


,Q, 


where chiral 


symmetry is restored in the nucleonic phase*. 

The separation of the chiral and the deconhnement transition persists also at 
non-zero T, where we checked two values T = 50 MeV and T = 100 MeV. However, 


It is important to stress this does not mean that the nucleons are massless, but that the parity 
partners are degenerate. 
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we do not expect the QMN model itself to be valid at such high T where gluons and 
lightest hadrons become important. Since in the QMN model we cannot calculate 
the Polyakov loop, the deconhnement transition in this case must be characterized 
in a different way. For these reasons, the £t of 14 to lattice Tc, performed here 
only in the very simple bQM model, carries systematic uncertainties. In general, 
the comparison to lattice QCD should be revised within a more complete approach 
before addressing the full T — fiB phase diagram. 

We must emphasize that the bag held should not be considered as an order pa¬ 
rameter in this model, since it is not connected to any of the fundamental QCD 
symmetries. In particular, it must not be considered as an order parameter for the 
deconhnement transition, even though it does play a key role in establishing it. The 
value of where the transition occurs is for some values of the a parameter ac¬ 
companied by a hnite jump in the total baryon density. In these cases, the density 
contrast between two coexisting phases can be used to characterize the deconhne¬ 
ment phase transition. We can draw an immediate analogy to the liquid-gas phase 
transition which does not have an order parameter related to a symmetry, but is 
also characterized by a jump in the density from the liquid to the gas phase. 

The approach to the deconhnement transition presented in this work is based 
on the phenomenological requirements and not rigorously grounded in QCD. Due 
to this, systematic uncertainties are expected, especially concerning the freedom in 
choice of the bag potential and the explicit form for the modihcation of the distri¬ 
bution functions. It would be very interesting to obtain a hrst-principle information 
about the hnite density distribution functions from Dyson-Schwinger studies, see 


e. g. 


73| . This microscopic input could then be used to construct more realistic 


ehective models of nuclear and quark matter. 
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VII. APPENDIX 


This appendix is devoted to solving a simplified version of the quark-meson- 
nucleon model, where we focus only on the role of the bag field b. The masses of 
nucleons and quarks are held fixed to 1000 MeV and 0 MeV, respectively. With 
all the other mean-fields except the b field discarded, the thermodynamic potential 
takes the following form 

d?p 


= 7x 

X=q,N± 


{2Tlf 


Tlog(l - nx) 
-hTlog(l - nx) 


(41) 


For the potential 14 we use the same Eq. fl25ll with the parameters Kb = 155 MeV, 
\b = 0.074. For the parameter a we take the value ab^ = 250 MeV. The partial 
pressures of quarks {pq) and nucleons {px) are calculated as 


Pq = -(H - Vbq) - n 


g > 


Pn — —^N+ — ^N- , 

where the pressure contribution arising from 14 is naturally assigned to quarks. 

From the minimization of the thermodynamic potential fjTTll we obtain 6 as a 
function of ps shown on Fig. [HI The b field follows its vacuum value up to ps = 1000 
MeV, then it begins to increase. The increase is due to the fact that nucleons 
favor a hnite value oi hb according to Eq. fl3^ . We expect the onset of quark 
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FIG. 8: The bag field as a function of in the simplified quark-meson-nucleon model, 
denoted with blue, dash-dotted line. For comparison we also plot the bag field from the 
full quark-meson-nucleon model for two different values of a calculated in Sec. |Vl 

degrees of freedom around /is — 3bo — 1700 MeV. An explicit calculation gives 
Hb = 1515 MeV. This point marks the sudden drop of the Afield as seen in Fig. |Hi 
The qualitatively similar characteristics are found in Sec. |V] in the complete model 
for abo = 300 MeV, also shown in Fig. |8l However, for higher abo = 440 MeV the 
transition becomes gradual and shifts to higher hb- 

On Fig. [9] we plot the partial pressures as a function of fiB- The nucleon pressure, 
given by the thick dashed red line, shows significant deviations from its ideal gas 
formula, shown by a thin red line. In particular, at fiB = 1515 MeV the nucleon 
pressure suddenly drops to zero value. The pressure arising from the quarks at first 
turns to negative values. The reason behind this is the increase of the b field in 
the region /i^ = 1000 — 1515 MeV. However, we stress that the total pressure of 
the system is always non-negative and continuous. After hb = 1515 MeV, nucleons 
disappear and the quark pressure becomes positive. Moreover, since quarks become 
the dominant degrees of freedom in the system, the value fiB = 1515 MeV marks 
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(MeV) 


FIG. 9: Normalized partial pressures of the nucleon and the quark sectors as a function 
of /is in the simplified quark-meson-nucleon model. The ideal nucleon pressure is plotted 
for comparison. 


the deconfinement transition point for this simplihed model. 
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